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Abstract. We prove comparison, uniqueness and existence results for viscos- 
ity solutions to a wide class of fully nonlinear second order partial differen- 
tial equations F(x, u, du, d 2 u) = defined on a finite-dimensional Riemannian 
manifold M. Finest results (with hypothesis that require the function F to be 
degenerate elliptic, that is nonincrcasing in the second order derivative vari- 
able, and uniformly continuous with respect to the variable x) are obtained 
under the assumption that M has nonnegative sectional curvature, while, if 
one additionally requires F to depend on a in a uniformly continuous man- 
ner, then comparison results are established with no restrictive assumptions 
on curvature. 



1. Introduction 

The theory of viscosity solutions to nonlinear PDEs on R" (and on infinite- 
dimensional Banach spaces) was introduced by M. G. Crandall and P. L. Lions in 
the 1980's. This theory quickly gained popularity and was enriched and expanded 
with numerous and important contributions from many mathematicians. We can- 
not mention all of the significant papers in the vast literature concerning viscosity 
solutions and Hamilton- Jacobi equations, so we will content ourselves with referring 
the reader to [4] and the references given therein. 

More recently there have been various approaches to extend the theory of vis- 
cosity solutions of first order Hamilton- Jacobi equations, and the corresponding 
nonsmooth calculus, to the setting of Riemannian manifolds. This is a natural 
thing to do, because many functions arising from geometrical problems, such as the 
distance function to a given set of a Riemannian manifold, are not differentiable. 
Also, many important nonlinear equations full of geometrical meaning, such as the 
eikonal equations, have no classical solutions, and their natural solutions, which in 
this case we think are the viscosity solutions, are not differentiable (if some readers 
disagree with our saying that viscosity solutions are the natural notion of solution 
for eikonal equations, they might change their mind if they have a look at the re- 
cent paper [5] , where the authors construct a 1-Lipschitz function u defined on the 
closed unit ball B of M. n ,n > 2, which is differentiable on the open ball B, and 
such that ||Vm(x)|| = 1 almost everywhere, but Vw(0) — 0; that is, the eikonal 
equation |Vit(x)|| = 1 in B, u — on dB, admits some exotic almost everywhere 
solutions which are everywhere differentiable and are very different from its unique 
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viscosity solution, namely the distance function to the boundary dB, which is not 
everywhere diffcrcntiablc but is much more natural from a geometric point of view). 

Mantegazza and Menucci [5] studied viscosity solutions to eikonal equations on 
Riemannian manifolds, in connection with regularity properties of the distance func- 
tion to a compact subset of the manifold. In [2] a theory of (first order) nonsmooth 
calculus for Riemannian manifolds (possibly of infinite dimension) was introduced 
and applied to show existence and uniqueness of viscosity solutions to Hamilton- 
Jacobi equations on such manifolds. Simultaneously, Ledyaev and Zhu [5] developed 
a (first order) nonsmooth calculus on finite-dimensional Riemannian manifolds and 
applied it to the study of Hamilton- Jacobi equations from a somewhat different 
approach, related to control theory and differential inclusions. 

The usefulness of nonsmooth analysis on Riemannian manifolds has been shown 
in [6] , where viscosity solutions are employed as a technical tool to prove important 
results in conformal geometry. 

However, to the best of our knowledge, no one has yet carried out a systematic 
study of second order viscosity subdifferentials and viscosity solutions to second 
order partial differential equations on Riemannian manifolds. 

In this paper we will initiate such a study by establishing comparison, uniqueness 
and existence of viscosity solutions to second order PDEs of the form 

F{x, u, du, d 2 u) = 

where u : M — > R and M is a finite-dimensional complete Riemannian manifold. 
We will study the Dirichlet problem with a simple boundary condition of the type 
u = f on dtt, where f2 is an open subset of M; and also the same equation, with 
no boundary conditions, on all of M. 

Let us briefly describe the results of this paper. We begin with the natural 
definition of second order subjet of a function u : M — > R, that is J 2 '~u(x) = 
{(d<p(x), d 2 ip(x)) : (f G C 2 (M, R), / — ip attains a local minimum at x}. This is a 
nice definition from a geometric point of view, but it would be complicated and 
uneconomic to develop a nonsmooth calculus exclusively based on this definition. 
It is more profitable to try to localize the definition through charts and then use the 
second order nonsmooth calculus on R™ to establish the corresponding results on M. 
However, second derivatives of composite functions are complicated, so not every 
chart serves this purpose, and we have to work only with the exponential chart. It 
is not difficult to see that ((, A) G J 2 '~u(x) if and only if (£, A) G J 2 '~ (uoexp x )(0). 

2 — 

When one turns to the limiting subjet J ' u(x) (defined as the set of limits 
of sequences (£ n ,A n ), where (£ n ,A n ) G J 2, ~u{x n ) and x n converges to x), things 
become less obvious but, with the help of a lemma which relates the second deriva- 
tives of a function <p : M — > R to those of the function ip = ip o exp^. (at points 

' 2 — 

near the origin in TM X ), one can still show that (C,A) G J ' u(x) if and only if 

(S,A) G J 2 '~( W oexpJ(0). 

By using this characterization we can extend Theorem 3.2 of [4] to the Riemann- 
ian setting. This kind of result can be regarded as a sophisticated nonsmooth fuzzy 
rule for the superdifferential of the sum of two functions, and is the key to the 
proof of all the comparison results in [4] and in this paper. The result essentially 
says that if u%, are two upper semicontinuous functions on M, ip is a C 2 smooth 
function on M x M, and we assume that u>(xi,X2) = u\(xi) + 1^2(2:2) — i f(xi,X2) 
attains a local maximum at (x±, £2), then, for each e > there exist bilinear forms 
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Bi G C 2 ((TM) Xi ,R), i = 1,2, such that 



() \ 2 + 

-(^(xi,X 2 ),B l G J 



for i = 1, k, and the block diagonal matrix with entries Bi satisfies 

\ + II V) I< ( * S ° 2 )<^ + ^ 2 . (*) 

where A = d 2 p(x 1 ,x 2 ) G C 2 s (TM Xl x TM i2 ,E). This is all done in Section 2 of the 
paper. 

In the case M = R n this result is usually applied with ip(x,y) = ^\\x — y\\ 2 , 
whose second order derivative is given by the matrix 



a 



I -I 
-I I 



When applied to vectors of the form (v,v) in M" x R n this derivative vanishes, 
which allows one to derive from (*) that B\ < Bi (as quadratic forms). This in 
turn provides a very general form of comparison result for viscosity solutions of the 
equation F(x, u, du, cPu) = in which the continuous function F is assumed to be 
degenerate elliptic (that is nonincrcasing in the variable d 2 u), strongly increasing 
in the variable u, and uniformly continuous with respect to x. 

The natural approach in the Ricmannian setting is then to consider <p(x,y) = 
^d(x,y) 2 , where d is the Riemannian distance in M. Two problems immediately 
arise. First, the function ip is not diffcrcntiable in general if the points x, y are not 
suitably close to each other. This is unimportant because, in the proof of the main 
comparison result, we only need ip to be C 2 smooth on a ball of small radius around 
a point xq which is the limit of two different sequences x a and y a , and we have to 
evaluate d 2 ip at the points (x a ,y a ). 

The second problem, however, is substantial. The second derivative of the func- 
tion ip is a quadratic form defined on TM X x TM y , and what we would like is that, 
when applied to a vector of the form (v, L xy v), where L xy is the parallel transport 
from TM X to TM y along the unique minimizing geodesic connecting x to y, this 
derivative is less than or equal to zero. This way condition (*) would imply that 
L X2Xl (B2) < B\, where L X2Xl (B 2 ) is the parallel transport of the quadratic form 
Bi from TM X2 to TM Xl along the unique minimizing geodesic connecting x 2 to x\, 
defined by 

(L X2Xl (B 2 )v,v) := (B 2 (L XlX2 v),L XlX2 v). 

And therefore we should be able to conclude that, if F is continuous, strongly in- 
creasing in the variable u, and degenerate elliptic (that is F(x, r, £, B) < F(x, r, £, A) 
whenever A < B), then a natural extension to X := {(x, r, A) : x G M, r e 1, ( g 
TM X ,A G C 2 S (TM X )} of the notion of uniform continuity of F(x,r, £, A) with re- 
spect to the variable x (namely, that 

\F(y, r, L xy (, L xy P) - F(x, r, (, P)\ -> uniformly as y ->• x, 

which wc abbreviate by saying that F is intrinsically uniformly continuous with 
respect to x) would be enough to show that comparison holds. 
However, as we will show in Section 3, one has that 

d 2 ip(x,y){v,L xy v) 2 < 
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for all v G TM X if and only if M has nonnegative sectional curvature. Therefore, 
with this choice of ip, one can get results as sharp as those in R™ only when one deals 
with manifolds of nonnegative curvature. Nevertheless, if the sectional curvature 
K of M is bounded below, say K > —K , then one can show that 

<P<p(x,y)(v, L xy vf < 2K d(x,y) 2 \\v\\ 2 

for all v G TM X , and by using this estimation it is possible to deduce that, if one 
additionally assumes that F satisfies a certain uniform continuity assumption with 
respect to the variables x and D 2 u of the kind "for every e > there exists S > 
such that d(x, y) < 5 and P — L yx Q < SI imply F(y, r, L xy (, Q) — F(x, r, (,P) < e 
for all C G TM*, P G C 2 S {TM X ), Q G C 2 (TM y ), r G R", then the comparison 
principle holds for the equation F — (either with the boundary condition u = 
on dQ, or with the assumption that M has no boundary and the functions u, v for 
which one seeks comparison are bounded). This is all shown in Sections 4 and 5. 

In Section 6 we see that Perron's method works perfectly well in the Ricmannian 
setting. For instance one can show existence of viscosity solutions to the equation 
u + G(x, du, d 2 u) — on compact manifolds under the same continuity assumptions 
on G as those that we require for comparison. 

In particular, we get the following: if M is a compact manifold and G is de- 
generate elliptic and uniformly continuous in the above sense, then there exists a 
unique viscosity solution of u + G(x, du, d 2 u) = on M . If one additionally assumes 
that M has nonnegative sectional curvature then the above uniform continuity as- 
sumption can be relaxed: it is enough to require that G is intrinsically uniformly 
continuous with respect to x, meaning that "for every e > there exists 5 > 
such that d(x,y) < 6 implies G(y, L xy (, L xy P) - G(x,(,P) < e for all ( G TM*, 

PeC 2 s {TM x y\ 

We end the paper by discussing the applicability of the above theory to some 
particular examples of equations. 

The notation we use is standard. M = (M, g) will always be a finite-dimensional 
Riemannian manifold. The letters X, Y, Z, V, W will stand for smooth vector fields 
on the Riemannian manifold M, and VyX will always denote the covariant de- 
rivative of X along Y. The Riemannian curvature of M will be denoted by R. 
Geodesies in M will be denoted by 7, a, and their velocity fields by 7', a' . If X is a 
vector field along 7 we will often denote X'(t) — %X(t) — V 7 i( t )X(t). Recall that 
X is said to be parallel along 7 if X'(t) — for all t. The Riemannian distance in 
M will always be denoted by d(x, y) (defined as the infimum of the lengths of all 
curves joining x to y in M). 

We will often identify (via the Riemannian metric) the tangent space of M at a 
point x, denoted by TM X , with the cotangent space at x, denoted by TM*. The 
space of bilinear forms on TM X (respectively symmetric bilinear forms) will be 
denoted by C 2 (TM X ) or £ 2 (TM X ,R) (rcsp. C 2 (TM X ) or C 2 (TM X , M)). Elements 
of C 2 (TM X ) will be denoted by the letters A,B,P,Q, and those of TM* by (,r], 
etc. Also, we will denote by T 2 , S (M) the tensor bundle of symmetric bilinear forms, 
that is 

T 2 , S (M)= |J C 2 S (TM X ,R), 

x£M 

andT 2;S (M) x =£2( TMx)R ). 
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We will make extensive use of the exponential mapping exp^, and the parallel 
translation along a geodesic 7 throughout the paper, and of Jacobi fields along 7 
only in Section 3. Recall that for every x 6 M there exists a mapping exp^,, defined 
on a neighborhood of in the tangent space TM X , and taking values in M, which 
is a local diffeomorphism and maps straight line segments passing through onto 
geodesic segments in M passing through x. The exponential mapping also induces 
a local diffeomorphism on the cotangent space TM*, via the identification given by 
the metric, that will be denoted by exp^ as well. 

On the other hand, for a minimizing geodesic 7 : [0, £\ — * M connecting x to y in 
M, and for a vector v G TM X there is a unique parallel vector field P along 7 such 
that P(0) = v, this is called the parallel translation of v along 7. The mapping 
TM X 3 » h P(f) g TM y is a linear isometry from TM X onto TM y which we will 
denote by L xy . Its inverse is of course L yx . This isometry naturally induces an 
isometry (which we will still denote by L xy ), TM* 3 ( h-> L xy ( G TM*, defined by 

(L xy (,v) y := ((,L yx v) x . 

Similarly, L xy induces an isometry C 2 (TM X , R) 3 A -> L xy (A) G £ 2 (TM y ,R) 
defined by 

(L xy (A)v,v) y := (A(L yx v),L yx v) x . 
By im{x) we will denote the injectivity radius of M at x, that is the supremum of 
the radius r of all balls B(0 x ,r) in TM X for which exp^, is a diffeomorphism from 
B(0 x ,r) onto B(x,r). Similarly, i(M) will denote the global injectivity radius of 
M, that is i(M) = mi{iM{x) : x <E M}. Recall that the function x 1— > Im{x) is 
continuous. In particular, if M is compact, we always have i(M) > 0. 

For Jacobi fields and any other unexplained terms of Riemannian geometry used 
in Section 3, we refer the reader to [3| [TO] . 

2. Second order viscosity subdifferentials on Riemannian manifolds 

Recall that the Hessian D 2 ip of a C 2 smooth function ip on a Riemannian man- 
ifold M is defined by 

D 2 ( p(X,Y) = (V x V<p,Y), 
where V</? is the gradient of ip and X, Y are vector fields on M (see [10] , page 
31). The Hessian is a symmetric tensor field of type (0, 2) and, for a point p G M, 
the value D 2 ip(X, Y)(p) only depends of / and the vectors X(p),Y(p) G TM p . 
So we can define the second derivative of tp at p as the symmetric bilinear form 
d 2 ip(p) : TM p x TM p -> R 

(v,w) ^ d 2 cp(p)(v,w) := D 2 ip(X,Y)(p), 

where X, Y are any vector fields such that X(p) — v, Y(p) = w. A useful way to 
compute d 2 ip(p)(v, v) is to take a geodesic 7 with 7'(0) = v and calculate 

d 2 

^V(7(*))|t=0, 

which equals d 2 ip(p)(v, v). We will often write d 2 ip(p)(v) 2 instead of d 2 (p(p)(v, v). 

Definition 2.1. Let M be a finite-dimensional Riemannian manifold, and / : M — > 
(—00, +00] a lower semicontinuous function. We define the second order subjet of 
/ at a point x G M by 

J '~f(x) = {(d(p(x) , d 2 (p(x)) : ip G C 2 (M, R), /— y attains a local minimum at x}. 
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If ((, A) £ J 2, ~/(x), we will say that £ is a first order subdifferential of /, and A is 
a second order subdifferential of / at x. 

Similarly, for an upper semicontinuous function g : M — > [—00, +00), we define 
the second order superjet of / at x by 

J ,+ f(x) = {(dip(x) , d 2 ip(x)) : <p £ C 2 (M,R), /— ip attains a local maximum at x}. 

Observe that J 2 ~f(x) and J 2 >+f(x) are subsets of TM* x C 2 (TM X ,R), where 
C 2 S (TM X ,R) = C 2 S {TM X ) stands for the symmetric bilinear forms on TM X . It is 
also clear that J 2 ~f(x) = — J 2 ' + (—f)(x), and that we obtain the same definitions 
if we replace the condition"^ £ C 2 (M, K)" with "93 is C 2 smooth on a neighborhood 
of x". 

By using the fact that a lower semicontinuous function / is bounded below on 
a neighborhood B of any point x with f(x) < 00, one can easily find a function 
ip £ C 2 (M, M) such that inf y6 as(/ — (p)(y) > f(x), hence f — f attains a local 
minimum at some point z £ B, and (d(p(z), d 2 ip(z)) £ J 2 '~f(z). This shows that 
the set 

{z £ M : J 2 '-f(z) ^ 0} 

is dense in the set {x £ M : /(a;) < 00}. A similar statement is true of upper 
semicontinuous functions. Therefore, when dealing with semicontinuous functions, 
one has lots of points where these subjets or superjets are nonempty, that is lots of 
points of second order sub- or super-differentiability. 

In the sequel M will always denote an n-dimensional Riemannian manifold. 
We next state and prove several results for subjets which also hold, with obvious 
modifications, for superjets. 

Proposition 2.2. Let f : M — ► (—00, +00] be a lower semicontinuous function. 
Let ( £ TM*,A £ C 2 (TM X ,M), x £ M. The following statements are equivalent: 

(1) (C,A)G J 2 --f(x). 

(2) /(exp») > f(x) + (C, v) x + \{Av, v) x + o(\\v\\ 2 ). 

Proof. (1) =^> (2): If ((, A) £ J 2 ~f(x), by definition there exists if £ C 2 (M,K) 
such that / — <p attains a local minimum at x and £ = dip(x), A — d 2 ip(x). We may 
obviously assume that ip(x) — f(x), so we have 

f(y) - v(y) > 

on a neighborhood of x. Let us consider the function h(v) = (p(exp x (vj) defined on 
a neighborhood of 0^ in TM X . We have that 

h(v) - h(0) + (dh(0),v) x + ^{d 2 h(0)v,v) x + o(\\v\\ 2 ). 

By taking y = cxp x (v) and combining this with the above inequality we get 

/(exp») > f(x) + (dh(0),v) x + ±(d 2 h(0)v,v) x +o(\\v\\ 2 ), 

so we only need to show that ( — dh(0) and A = d 2 h(0). To see this, let us fix 
v £ TM X and consider the geodesic -f(t) — exp x (tv) and the function t ^ ip(j(t)) = 
h(tv). We have that 
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and 

^h{tv) = {d\{ 1 {t)) 1 '{t), 1 '{t)). 
In particular, for t = 0, we get 

dh(0)(v) = jh(tv)\ t=0 = (d<p(x),v) = (C,v), 

that is dh(0) = £; and also 

d 2 

(d 2 h{0)v,v) = —h(tv)\t= = (d 2 ip(x)v,v) = (Av,v), 
that is A = d 2 h(0). 

(2) => (1): Define F(v) = f(exp x (v)) for v in a neighborhood of CL, G TM X . We 
have that 

F(v) > F(0) + (C, v) x + \{Av, v) x + o(\\v\\ 2 ). 

The result we want to prove is known to be true in the case when M — R™, so 
there exists ip : TM X — ► R such that F — ip attains a minimum at and dip(0) = (, 
d V'(O) = ^4- Since minima are preserved by composition with diffemorphisms, the 
function ip := V ° exp^ x , defined on an open neighborhood of a; G M, has the 
property that f — tp = (F — ip) o exp^T 1 attains a local minimum at x = exp~ 1 (0). 
Moreover, according to (1) =^> (2) above, we have that 

dkp{x) = #(0), and d 2 (p(x) = d 2 iP{0), 

so we get dip(x) = ( and d 2 (p(x) = A. Finally, by using smooth partitions of unity 
we can extend ip from an open neighborhood of x to all of M. □ 

Corollary 2.3. Let f : M — ► (— oo, +oo] be a lower semicontinuous function, and 
consider ( G TM* , A G C 2 S (TM X ,R), x G M . Then 

(C) A) G J 2 '~f(x) (C,A)G J 2 --(/oexp x )(0 x ). 

Making use of the above characterization, one can easily extend many known 
properties of the sets J 2 '~ f(x) and J 2 ' + f(x) from the Euclidean to the Riemannian 
setting. For instance, one can immediately see that J 2, ~f(x) and J 2 ' + f(x) are 
convex subsets of TM* x C 2 (TM X ). They arc not necessarily closed, but if one 
fixes a C G TM* then the set {A : (£,A) G J 2 '~f(x)} is closed. A useful property 
that also extends from Euclidean to Riemannian is the following: if ip is C 2 smooth 
on a neighborhood of x then 

J 2 -(/-^)(x)={(C-^(x), A-d 2 ^{x)) : (C,A)eJ 2 --f(x)}. 

One can also see that / is twice differcntiablc at a point x G M (in the sense 
that for some (unique) C S TM*,A G £ 2 S (TM X ,R) we have that f(exp x (v)) = 
f(x) + (C,v) + \{Av,v) + o(\\v\\ 2 ) asu^O) if and only if J 2 ~f{x) D J 2 ' + .f{x) is 
nonempty (in which case J 2 '~f(x) n J 2 ' + f(x) = {((,A)}). 

Next we have to define the closures of these set-valued mappings. Let us first re- 
call that a sequence (A n ) with A n G C 2 (TM Xn ) is said to converge to A G C 2 (TM X ) 
provided x n converges to x in M and for every vector field V defined on an open 
neighborhood of x we have that (A n V(x n ), V(x n )) converges to (AV(x),V(x)). 
Since we have (AV, W) = \ ({A(V + W), V + W) - (AV, V) - {AW, W)), it is clear 
that this is equivalent to saying that (A n V(x n ), W{x n )) converges to (AV(x), W{x)) 
for all vector fields V, W on a neighborhood of x in M. 
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Similarly, a sequence (£„) with £„ € TM*^ converges to £ provided that x n — > x 
and (£ n , V(x n )) —> ((, V(x)) for every vector field V defined on an open neighbor- 
hood of x. 

Remark 2.4. It is not difficult to see that, if M = R", then A n (respectively £„) 
converges to A (resp. C) m the above sense if and only if \\A n — A\\ — > (resp. 
llCn - Cll 0) m £ 2 (R n ,R) (resp. in R"). 

It is also worth noting that \\A n - A\\ (resp. ||£„ - Cll -» 0) in £ 2 (R",R) 
(resp. in R") if and only if (A n v, v) — > (Av,v) (resp. {( n ,v) — > ((,v)) for every 
v G R", that is pointwise convergence is equivalent to uniform convergence on 
bounded sets, as far as linear or bilinear maps on 1" are concerned. 

Definition 2.5. Let / be a lower semicontinuous function defined on a Riemannian 
manifold M, and x 6 M. We define 

J 2 '~/(.t) = {(C,A) e TM* x £ S (TM X ) : 3x„ e M,3(C„, A„) e J 2 - f(x n ) 

s.t. (x n , f(x n ),(n,A n ) -> (a;, /(x), C, A)}, 

2.+ 

and for an upper semicontinuous function g on M we define J ' (/(a;) in an obvious 
similar way. 

Remark 2.6. According to Remark 12.41 we have that, in the case M = R™, the 

2 — 2 + 

sets J ' g{x) and J ' g(x) coincide with the subjets and superjets defined in [4]. 

2 — 

In order to establish the analogue of Corollary 12.31 for the closure J ' g{x), we 
will use the following fact. 

Lemma 2.7. Let tp : M — > R be a C 2 smooth function, and define ip = (poexp x on a 
neighborhood of a point € TM X . Let V be a vector field defined on a neighborhood 
of in TM X , and consider the vector field defined by V(y) = dexp x (w y )(V(w y )) 
on a neighborhood of x in M, where w y :— exp~ 1 (j/), and let 

a y (t) = exp x (w y +tV(w y )). 

Then we have that 

D 2 iP(V, V)(w y ) = D 2 V {V, V)(y) + (V<p(y), <(0)). 
Observe that c"(0) = so, when y — x, we obtain 

d 2 i/j(0){v,v) = d 2 ip{x){v,v) 

for every v € TM X . 

Proof. Fix y near x. We have that 

j t iP{w y +tV{w v )) = j t v{e y {t)) = (Vv(o- y (t)),o-' y (t)}, 

and 

^ip(w y +tv(w y )) = ^KM) = (v CT , (t) v^K(i)),<W) + (v^K(i)),<W). 

Note that o-' y (0) — V(y), hence by taking t — we get the equality in the statement. 
Observe that when y — x the curve a x is a geodesic, so a"(0) =0. □ 
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Proposition 2.8. Let f : M — ► (— oo, +00] be a lower semicontinuous function, 
and consider ( G TM*,A G C 2 (TM X , R), x G M. Then 

((,A)eJ 2 '-f(x) ^ (C,A)€ J 2 '"(/oex Px )(0 x ). 

Proof. (=>) : If (C-A) £ </ 2 ' /(#) there exist x n — > ir and (Cn, Ai) € J 2 '~f(x n ) so 
that Cn - * C) ^« — > -4, /(^n) — * /(^O- Take G C 2 (M) such that / — ip n attains 
a minimum at x n and Cn = dip n (x n ), A n = d ip n (x n ). Define ip n = ip n o exp^, on a 
neighborhood of in TM XJ and v n = exp~ It is clear that foexp x —ip n attains 

a minimum at v n . We then have that {dip n {v n ) 1 d 2 '^ n {v n )) g J 2, ~(f o exp x )(v n ), 
and since t>„ — * and /oexp 2 ,(i; n ) — > /(a;), we only have to show that dip n (v n ) — > £ 
and d 2 ipn{v n ) — » -A- 

Take a vector field V" on TM X , and define a corresponding vector field V on a 
neighborhood of x in M by 

V(y) = dexp x (w y ){V{w y )), 
where w y = exp~ 1 (y). We have that 

(dtp n (v n ),V(v n )) = (dtp n (x n ) o dcxp x (v n ),V(v n )) = (dip n (x n ),V(x n )), 
so we get 

(#„K),FK)) - (C n ,v(x n )) -> (C, = (C,no)), 

which shows dip n (v n ) — ► C- O n the other hand, according to the preceding Lemma, 
we also have that 

dfynK)(V(«„), K(« n )) = A n (V(x n ), V{x n )) + (Cn, <,(0)), 

where <Tj,(t) = exp^tt;^ + 

Notice that the mapping y \— ► (7^(0) defines a smooth vector field on a neighbor- 
hood of x in M (and in particular <J Xn (0) — > c"(0) = as n — > 00). Since A„ — > A, 
C„ — ► C, we get, by taking limits as n — > 00 in the above equality, that 

dVn(«n)(^K), V(w„)) -» A(V(x), V(x)) + = A(V(0), V(0)), 
which proves that d 2 tp n (v n ) — ► A. 

(-4=) If (C,A) g o exp x )(0) there exist u n -> and {C, n ,A n ) g J 2 '~(/ o 

ex Pj(«n) so tnat Cn C> A» ~* A, f(x n ) -> /(a;), where ar n = exp x (w n ). J^ake 
VVi g C 2 (TM X ) such that / o exp^ — ^„ attains a minimum at v n and Cn = 
dip n (vn), A n = d 2 ip(v n ). Define tp n = tp n o cxp^ 1 on a neighborhood of x in M. 
Then / — ip n attains a minimum at x n , so (dip n (x n ), d 2 ip n (x n )) G J 2 '~ f(x n ), and 
we only have to show that dip n (x n ) — > C and d 2 ip n (x n ) — > A. Take a vector field 
V on a neighborhood of x in M, and define a corresponding vector field V on a 
neighborhood of in TM X by 

V(w v ) = dexp- 1 (y)(V(y)), 

where w y = cxp x 1 (y). Now we have that 

(dip n {v n ), V(v n )) = (dip n (x n ),V(x n )), 
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from which we deduce that dip n (x n ) — » Ci an( i also, by using this fact and the 
preceding Lemma, 

d 2 (fi n (V(x n ), V(x n )) = 

A n (V(v n ),V(v n )) - (&p n (x n ),oZ n (0)) -» 

- (4v(i), v(x)) - <c,4'(o)) = -o, 

concluding the proof. □ 

Remark 2.9. One can see, as in the case of J 2 '~ f(x), that if ip is C 2 smooth on 
a neighborhood of x then 

J 2 '-(/-^)(z) = {(C A-dfyOc)) : ((,4)eJ 2 '7(4 

The following result is the Ricmannian version of Theorem 3.2 in [4] and, as in 
that paper, will be the key to the proofs of comparison and uniqueness results for 
viscosity solutions of second order PDEs on Ricmannian manifolds. 

Theorem 2.10. Let M±, M}. be Riemannian manifolds, and flj C Mi open 
subsets. Define Q = Oi x . . . x fife C Mi x . . . x Mfc = M . Let Ui be upper 
semicontinuous functions on £li, i = 1, ...,k; let ip be a C 2 smooth function on 
and set 

lu(x) = ui(xi) + . . . + u k (xk) 

for x = (x±, ...,Xk) £ ^- Assume that x — (x±, . . . , &k) is a local maximum ofuj — ip. 
Then, for each e > there exist bilinear forms Bi € C 2 ((TMi) Xi ,M), i — 1, k, 
such that 

^-<p(x),B t ^J S J 2 ' + u i (x l ) 

for i = 1, k, and the block diagonal matrix with entries Bi satisfies 

Bx ... 



\A\\ \I< I : : \<A + eA 2 , 



. . . B k 



where A = d 2 ip(x) € C 2 {TM X ,} 



Recall that, for ( e TM*, A e C(TM X x TM X ,R), the norms ||C|U an d \\A\\ X 
are defined by 

IICIU = sup{(C,u)x : v 6 TM„, ||u|| x < 1}, 

and 

= sup{\(Av,v) x \ : v E TM X , \\v\\ x < 1} = sup{|A| : A is an eigenvalue of ^4}. 

Proof. The result is proved in in the case when all the manifolds Mi are Euclidean 
spaces, and we are going to reduce the problem to this situation. By taking smaller 
neighborhoods of the Xi if necessary, we can assume that the fii are diffeomorphic 
images of balls by the exponential mappings exjp x . : B(0, r^) — ► fij = B(ii, r^), and 
that exp s maps diffeomorphically a ball in TM X onto a ball containing fl. The 
exponential map exp a from this ball in TM X = (TM\) Xl x . . . x (TMk)x k into M 
is given by 

exp £ («i, ...,v k ) = (exp^Oi), ...,exp £fc (i; fe )) . 
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Now define functions on open subsets of Euclidean spaces by ui(v) = cj(exp £ (u)) 
and Ui{vi) = Ui(exp x .(vi)). We have that u>(vi, ...,Vk) — Si(ui) + ... + Uk(vk), and 
Ox — (Oij , Xk ) is a local maximum of ui — ip, where %p = ip o exp £ . 

Then, by the known result for Euclidean spaces, for each e > there exist 
bilinear forms B. t e C 2 s ((TMi) Xi ,R), i = 1, k, such that 

JLi>(0 & ),b?) eJ 2 ' + ^(0 £i ) 



9«i 

for i = 1, k, and the block diagonal matrix with entries Bi satisfies 

\ + \\AU I< : : I .1+ 

V ... B k 

where A = d 2 ip(O x ) G Z^fjTM^R). According to Proposition 12 . 81 we have that 

^(O i ),B,jeJ u 5,(0 ii ) <s=^ ^A^),^ g J 2,+ ^(£ 4 ), 

so we are done if we only see that 

g 
jr- iPi^x) = -^-<f(x) and d 2 ifj(O x ) = d 2 ip{x). 



But this is a consequence of Lemma 12.71 □ 



Now we extend the notion of viscosity solution to a Hamilton- Jacobi equation 
on a Ricmannian manifold. In the sequel we will denote 

X := {(x, r, C, A) : x G M, r G R, C G ™», A 6 ^(TM X )} 

Definition 2.11 (Viscosity solution). Let M be a Riemannian manifold, and F : 
X — > R. We say that an upper semicontinuous function u : M — > R is a viscosity 
subsolution of the equation F — provided that 

F(x,u(x),t,A) <0 

for all cc G M and (C, -A) G J 2 ' + u(a;). Similarly, a viscosity supersolution of F = 
on M is a lower semicontinuous function u : M — > R such that 

F(x,u(x),£,A) > 

for every x G M and (C A) G J 2 '~u(a;). If u is both a viscosity subsolution and a 
viscosity supersolution of F = 0, we say that u is a viscosity solution of F = on 
M. 

Remark 2.12. If u is a solution of F < and F is continuous on ,Y then 

2,+ 

F(x, u(a;), £, A) < for every ((, A) G J ' A similar observation applies 

to solutions of F > and solutions of F = 0. 

Definition 2.13 (Degenerate ellipticity). We will say that a function F : X — > R 
is degenerate elliptic provided that 

A<S F(x,r,C,B) <F(x,r,C,A) 

for all a; G M, r G R, C € TM X , A,Be C 2 (TM X ). 
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Example 2.14. If we canonically identify the space of symmetric bilinear forms 
on TM X with the space of self-adjoint linear mappings from TM X into TM X , we 
have that 

Lyx{Q) = L X yQL xy . 

Hence 

tra,ce(L yx Q) = trace(Q), and det + (L yx (Q)) = det+(Q) 

(where det + A is defined as the product of the nonnegative eigenvalues of A), and 
it is immediately seen that the functions G(r,(,A) = —det + (A) and H(r, (,A) — 
— trace(A) are degenerate elliptic and, moreover, are invariant by parallel transla- 
tion, in the sense that 

G(r,(,A)=G(r,L xy (,L xy A) 

for all r E K, £ £ TM X ,A 6 C 2 (TM X ). The same can be said of all nonincreasing, 
symmetric functions of the eigenvalues of A. Thus one may combine such functions 
to construct many interesting examples of equations to which our results apply, as 
we will see later on. 

Remark 2.15. If the function F is degenerate elliptic, then every classical solution 
of F = is a viscosity solution of F = 0, as is immediately seen. However this may 
be not true if F is not degenerate elliptic; for instance when M = R the function 
u(x) — x 2 — 2 is a classical solution of u"(x) + u(x) — x 2 = but is not a viscosity 
solution. 

In order that the theory of viscosity solutions applies to an equation F — 0, the 
following condition is usually required. 

Definition 2.16 (Properness). We will say that a function F : X — > K., (x, r, £, A) i— » 

F(x, r, C A), is proper provided 

(i) F is degenerate elliptic, and 

(ii) F is nondecreasing in the variable r. 



3. A KEY PROPERTY OF THE HESSIAN OF THE FUNCTION (x, y) H-> d(x, lj) 2 

When trying to establish comparison results for viscosity solutions of second 
order PDEs on a Riemannian manifold M we will need to know how the Hessian 
of the function ip : M x M — > E, 

<p(x,y) = d(x,y) 2 

behaves. More precisely we will need to know on which manifolds M one has that 

d 2 ^(x,y){v,L xy v) 2 <0 (|t) 

for all v £ TM X) with x,y E M close enough to each other so that d(x,y) < 
mm{i M (x), iu (y)}- 

Let us calculate this derivative. We have that 

^(x,y) = 2d(x,y)^(x,y) = -2exp- 1 ( 2/ ). (1) 

The second equality can be checked, for instance, by using the first variation formula 
of the arc length (see [101 p. 90]). Indeed, if a(t, s) is a variation through geodesies 
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of a minimizing geodesic j(t) with y = 7(0) and x = where I = d(x, y), and if 
L(s) denotes the length of the geodesic 1 1— > a(t, s), then 



(V,T)f - / (V,V T T)dt 



((V(£),T(l)) (V(0),T(0)» 



where T = da/dt (so VtT = 0) and V = da/ds. Taking an a such that V is the 
Jacobi field along 7 satisfying F(0) = 0, V(£) = v, we get 



^(x,y)(v) = ^-L(s)\ s = = -^-exp^y)). 



Similarly, we have 



-^{x,y) = 2d{x,y) — {x,y) = -2cxp- 1 (x). (2) 



Observe that 



dip. . T .dp , 9d. . T ,dd . , . 

— + L xy{-Q^ (x,y)) = = Q^{x,y) + L xy { — {x,y)). (3) 

By differentiating again in (1) and (2), we get 

00c, y)(«) 2 - 2 (^(x, y)( V )) + 2d(x, y)(i0 2 , 

<9 2 (z> , .. . „<9d 9t! <9 2 <i . .. . 

10) = 2— (x,y)(v) — (x,y)(w) + 2d(a;, y) — — -(a;, y)(v, w) 



dxdy ' dx ' dy ' ' ctedy 

|^(x,y)H 2 =2 ^(x,y)H^) + 2d(x,y)|-^(x,y)H 2 , 

so, if we take w = L xy v and we sum the two first equations, and then we use (3), 
we get that 

-^(x,y){v) 2 +-^-(x,y){v,L xy v) = 2d{x,y) 



d 2 d, w , 2 d 2 d 



dx^ y){v) + dx^-y {x ^ LxyV \ 

and we get a similar equation by changing x for y. By summing these two equations 
we get 

d 2 ip(x,y)(v,L xy v) 2 = 2d(x,y)d 2 (d)(x,y)(v,L xy v) 2 , 
so it is clear that condition (jj) holds if and only if 

d 2 (d)(x,y)(v,L xy v) 2 <0 (b) 

for all v G TM X . 

Another way to write conditions (jj) or (b) is 

d 2 
dt 2 

where a x and a y are geodesies with a x (0) = x, cr y (0) — y, <j' x (0) = v and a' y (0) = 
L xy v. The function t ^ h{t) := d(a x (t),a y (t)) measures the distance between 
the geodesies a x and a y (which have the same velocity and are parallel at t — 0) 
evaluated at a point moving along any of these geodesies. 

We are going to show that the second derivative h"(0) is negative (that is, con- 
dition (t]) holds) if and only if M has positive sectional curvature. 

In particular, by combining this fact with Equation (3) (which tells us that 
h'(0) = 0), we see that the function h(t) attains a local maximum at t = if 



{d(<T x (t),a y (t)))\ t =o < 0, 
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and only if M has positive sectional curvature. This corresponds to the intuitive 
notion that two geodesies that are parallel at their starting points will get closer 
if the sectional curvature is positive, while they will spread apart if the sectional 
curvature is negative. 

Proposition 3.1. Condition (ft) (equivalently (b), or (§)) holds for a Riemannian 
manifold M if and only if M has nonnegative sectional curvature. In fact one has, 
for the function <p(x, y) = d(x, y) 2 on M x M , that: 

(1) If M has nonnegative sectional curvature then 

d 2 ip(x,y)(v,L xy v) 2 < 

for all v G TM X , with x,y G M close enough to each other so that d(x, y) < 
mm{i M (x),i M (y)}- 

(2) If M has nonpositive sectional curvature then 

d 2 <f{x,y)(v,L xy v) 2 > 

for all v G TM X , x,y G M such that d(x,y) < mm{iM{x), (?/)}• 

This fact must be known to the specialists in Riemannian geometry, but we have 
been unable to find a reference for part (1), so we provide a proof. Let us begin by 
reviewing some standard facts about the second variation of the arc length and the 
energy functionals. 

Take two points Xo,yo G M with d(xo,yo) < minjzM^o), *m(?/o)}j an d let 7 be 
the unique minimizing geodesic, parameterized by arc- length, connecting x to y . 
Denote I — d(xo,yo), the length of 7. Consider a(t,s), a smooth variation of 7, 
that is a smooth mapping a : [0,£] x [— e, e] — ► M such that a(t, 0) = "f(t) for all 
t G [0,1]- Consider the length and the energy functionals, defined by 

L(s) = L{a s )= f\\a' s {t)\\dt 
Jo 

and 

E(s) = E(a s )= [ \\a' s (t)\\ 2 dt, 
Jo 

where a s is the variation curve defined by a s (t) = a(t,s) for every t G [0,i\. 
According to the Cauchy-Schwarz inequality (applied to the functions / = 1 and 
g(t) = ||ag(i)|| on the interval [0,£]) we have that 

L(s) 2 < lE(s), 

with equality if and only if ||o4(t)|| is constant. Therefore, in the case when a s is 
a geodesic for each s (that is a is a variation of 7 through geodesies) we have that 

L(s) 2 = £E(s) 

for every s G [— e, e]. 

Now take a vector v G TM Xo , set w — L XoVo v, and consider the geodesies a X(i ,a yo 
defined by 

a X0 (s) = exp xa (sv), a yo (s) = exp yo (sw). 

We want to calculate 

d 2 

d 2 ip(x„,y )(v,w) 2 = —(p(a xo (s),a yo (s))\ s=Q , 
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where ip(x,y) = d(x,y) 2 . To this end let us denote by a s : [0,£] — > M the unique 
minimizing geodesic joining the point a Xo (s) to the point <J Vo {s) (now, for s ^ 0, 
a s is not necessarily parameterized by arc- length), and let us define a : [0,£] x 
[— e, e] — > M by a(t, s) — a s (t). Then a is a smooth variation through geodesies of 
j(t) = a(t, 0) and, according to the above discussion, we have 

ip(<7 Xo (s),a yo (s)) =L(s) 2 =£E(s), 

and therefore 

d 2 y{x ,y ){v,wY =IE"(0). (4) 

If we denote X(i) = da(t, 0)/ds, the variational field of a, then the formula for the 
second variation of energy (see O p. 197]) tells us that 

\e"(0) = -J o (X,X''+i?(^X) 7 0dt+(X(i) ) X'( i ))|g + (^(t,O) )T '(f))|^ o , 

(5) 

or equivalently 

\E"{Q) = f ((X>, X') {R{i, X)j>, X)) dt + (~(t, 0), j'(t)) \tzi (6) 

where we denote X' = V 7 '( t )X, and X" = V 7 /( t )X'. 

Note that, since the variation field of a variation through geodesies is always a 
Jacobi field, and since the points xq and yo are not conjugate, the field X is in fact 
the unique Jacobi field along 7 satisfying that X(0) = v, X(£) — w, that is X is 
the unique vector field along 7 satisfying 

X"(t) + R(-y'(t),X(t))i'(t) = 0, and X(0) = v, X{t) = w, 

where R is the curvature of M. On the other hand, since the curves s — > a(0, s) — 
a Xo (s) and s — ► a(^, s) = cr ya (s) are geodesies, we have that 

(^M),y»>iS = o. 

These observations allow us to simplify (5) and (6) by dropping the terms that 
vanish, thus obtaining that 

±E"(0) = {X(£),X'(e))-{X(0),X'(0)) (7) 

and also 

l -E"{Q) = J ((X', X') <i?( 7 ', XW, X)) dt. (8) 

Recall that the right-hand side of (8) is called the index form and is denoted by 
I(X,X). 

By combining (4), (7) and (8) we get 

d 2 ^(x ,y )(v,w) 2 = 2£((X(£),X\£)) - (X(Q),X'(0))) = 

= 21 f ((x\x') - (R(y,x)y,x))dt. 

Jo 

Therefore condition (fl) holds if and only if, for every Jacobi field X along 7 with 
X(0) = v,X(£) = w = L Xoyo v, one has that 

{X(£),X'(£))-{X(0),X'(0))<0 (0) 
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or, equivalently, 

/ i(X', X') - (Rh', X)i, X)) dt < 
Jo 

for the same Jacobi fields. 

Proof of Proposition I3.lt The proof of (2) is immediate and is well referenced 
(see for instance [7J Theorem IX.4.3]): if M has nonpositive sectional curvature 
then we have (R^ 1 , X)j' , X) < 0, hence, according to the above formulas, 

dV(*o, Vo)(v, w) 2 = 2£ / ({X 1 , X 1 ) - {R{i, X)i, X)) dt>2£ (X 1 , X')dt > 0, 



which proves (2). Note that in this case we do not use that w = L Xoyo v 7 so this 
holds for all v,w. 

Our proof of (1) uses the following Lemma, which is a restatement of Corollary 
10 in Chapter 8 of QI]. 

Lemma 3.2. Let 7 : [0,1] — » M be a geodesic without conjugate points, X a Jacobi 
field along 7, and Z a piecewise smooth vector field along 7 such that X(0) = Z(0) 
andX(£) = Z(£). Then 

I(X,X)<I(Z, Z), 
and equality holds only when Z = X . 

That is, among all vector fields along 7 with the same boundary conditions, the 
unique Jacobi field along 7 determined by those conditions minimizes the index 
form. Recall that 

I(Z, Z) = I ((Z\ Z') - (R( 7 ', Z)i, Z))dt, 



but this number is not equal to {Z(t),Z'{t)) - (Z(0),Z'(0)) unless Z is a Jacobi 
field. 

Let X be the unique Jacobi field with X(0) — v,X(£) — w = L Xoyo (v). Define 
Z = P(t), where P(t) is the parallel translation along 7 with P(0) — v (hence 
P(£) = w). The field Z is not necessarily a Jacobi field, but it has the considerable 
advantage that Z'(t) = for all t, so we have that 

I(Z,Z)= [ ((Z',Z')-{R(i,Z)i,Z))dt = - I (R{i,Z)i,Z)dt<Q 



because M has nonnegative sectional curvature. We then deduce from the above 
Lemma that 

(X(£), X'{£)) - (X(0), X'(0)) - I(X, X) < I(Z, Z) < 0, 
which, according to the above remarks (see (<^)), concludes the proof. □ 

Even though we will not have d 2 ip(x, y)(v, L xy v) 2 < when M has negative 
curvature, we can estimate this quantity and show that it is bounded by a term 
of the order of d(x,y) 2 , provided that the curvature is bounded below. This will 
also be used in the next section to deduce a comparison result which holds for all 
Riemannian manifolds (assuming that F is uniformly continuous). 
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Proposition 3.3. Let M be a Riemannian manifold. Consider the function tp(x, y) = 
d{x,y) 2 , defined on MxM. Assume that the sectional curvature K of M is bounded 
below, say K > —Kq. Then 

d 2 tp(x,y)(v,L xv v) 2 < 2K d(x,y) 2 \\v\\ 2 

for all v £ TM X and x,y G M with d(x,y) < min{iM(x), iu (y)}- 

Note that for Kq — we recover part (1) of Proposition 13. II 



Proof. Let X, Z be as in the proof of (1) of the preceding Proposition. With the 
same notations, we have that 

d 2 <p(x,y)(v,L xy v) 2 = 2£I(X,X) < 2£I(Z,Z) = 

-21 I (R(~f',Z)~f',Z)dt <2£ K \j'(t) AZ(t)\ 2 dt < 
Jo Jo 



2t Ko||7'WI| 2 ||^WI| 2 = 2^ / if |M| 2 cft = 2^ |M| 2 , 
Jo Jo 

which proves the result. □ 



4. Comparison results for the Dirichlet problem 

In this section and throughout the rest of the paper we will often abbreviate 
saying that u is an upper semicontinuous function on a set fl by writing u G 
USC(Ci). Similarly, LSC(Cl) will stand for the set of lower semicontinuous functions 

on n. 

The following lemma will be used in the proof of the main comparison result for 
the Dirichlet problem 

F(x,u(x),du(x),d 2 u(x)) = on ft; u = f on 3ft. {DP) 

Lemma 4.1. Let ft be a subset of a Riemannian manifold M, u 6 USC(Q), 
v e LSC(H) and 

m a := sup(w(x) - v(y) - -d(x,y) 2 ) 

fixf! ^ 

for a > 0. Let m a < oo for large a and (x a , y a ) be such that 

lim (m a - (u(x a ) - v(y a ) - ^rd{x a ,y a ) 2 )) = 0. 

a — *oo Z 

Then we have: 

(1) lim Q ^oo ad(x a , y a ) 2 — 0, and 

(2) liniQ^oo m a = u(x) — v(x) — sup x ^a(u(x) — v(x)) whenever x G ft is a 
limit point of x a as a — > oo. 

Proof. The result is proved in [U Lemma 3.1] in the case M = M. n , and the same 
proof clearly works in the generality of the statement (in fact this holds in any 
metric space). □ 

Now we can prove the main comparison result for the Dirichlet problem. 

Theorem 4.2. Let ft be a bounded open subset of a complete finite- dimensional 
Riemannian manifold M , and F : X — > M. be proper, continuous, and satisfy: 
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(1) there exists 7 > such that 

7 (r - s) < F(x,r,C,Q) - F(x,s,C,Q) 

for r > s; and 

(2) there exists a function lu : [0, 00] — > [0, 00] with lim f ^ + ^(t) = and such 
that 

F(y,r,aexpy 1 (x),Q) - F(x,r,-aexp~ x (y),P) < uj(ad(x,y) 2 + d(x,y)) 
for all x,y G n, r G R, P G T 2 , a {M) x , Q g T 2 , s (M) y with 

-{i + M ){o ?)^(o " Q )<A a+ s a Al, (*) 

where A a is the second derivative of the function tp a (x,y) = ^d(x,y) 2 at 
the point (x, y) g M x M , 

1 

£a ~ 2(1 + \\A a \\y 

and the points x, y are assumed to be close enough to each other so that 

d(x,y) < mhi{i M {x),i M (y)}- 
Let u g USC(fl) be a subsolution and v G LSC(fl) a supersolution of F = on Q, 
and u < v on d£l. 

Then u < v holds on all of O. 

In particular the Dirichlet problem (DP) has at most one viscosity solution. 

Proof. Assume to the contrary that there exists z G f2 with u(z) > v(z). By 
compactness of Q and upper semicontinuity of u — v, and according to Lemma |4~T| 
there exist x a , y a so that, with the the notation of Lemma |4. 11 

u(x a ) - v(y a ) - ^d(x a ,y a ) 2 = m a > 6 := u(z) - v(z) > 0, (3) 

and 

ad(x a , y a ) 2 — > as a 00. (4) 
Again by compactness of fl we can assume that a subsequence of (x a , y a ), which we 
will still denote (x a , y a ) (and suppose a g N), converges to a point (xq, yo) g Cl x Q. 
By Lemma 14. II we have that xq = yo and 

S < lim m a — u(xq) — v(xq) — sup(u(x) — v(x)), 

and in view of the condition «<uon dtt we have that xq G f2, and x a , y a G f2 for 
large a. 

Fix ro > and Rq > such that, for every x G B(xo,ro), exp^ is a dif- 
feomorphism from B(0,Rq) C TM x onto B(x,Rq) D B(xo,r ) (see [3j Theo- 
rem 3.7 of Chapter 3]). Then, for every x,y G B(xo,ro) we have that d(x,y) < 
min{iM(x) 7 «j\/(y)}, the vectors exp~ 1 (y) G TM X = TM* and exp~ 1 (a;) G TM y = 
TM* are well defined, and the function tp(x, y) = d(x, y) 2 is C 2 smooth on B(xq, ro) x 
B(xo,ro) G M x A/. Taking a subsequence if necessary, we can assume that 
x a ,ya g B(xo,ro) for all a. 

Now, for each a, we can apply Theorem 1 2 . 101 with Qi = Q 2 = B(xQ,ro), u\ = u, 
u 2 = -v, ip(x, y) = <p a (x, y) := fd(x, y) 2 , and for 

1 

£ - £ «'- 2(1 + 11^(^,^)11)- 
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Since (x a ,y a ) is a local maximum of the function (x,y) i— > u(x) — v(y) — tp(x,y), 
we obtain bilinear forms P E ^((TM Io ,I), and Q e £*((TAf"j, a ,R) such that 

-^ip(x a ,y a ),P^j EJ 2 ' + u(x a ), 
-J^¥>(£a,2/a),Qj 6 J 2: v(j/o) 

2 — 2.+ 

(recall that J ' f(y Q ) = —J ' (— v)(2/c*)), and 

j- + \\A a \\\l< ( ? ° Q \ <A a + e a Al 

where A a = d 2 (p(x a ,y a ) E C? s (TMr x ^ y \ , R) , so we get that condition (*) holds for 
.x = x a , y = y a . Therefore, according to condition (2), we have that 

F(y a ,r,aexp~*(xa,), Q) - F(x a ,r, -aexp~^(y a ), P) < uj(ad(x a ,y a ) 2 + d(x a ,y a )) 
On the other hand, from equation (1) in the preceding section we have that 
-^(p(x a ,y a ) = -aexp~^(y a ), and - -^-ip(x a ,y a ) = aexp~*(x a ), 

2 -4- 2 

hence (— aexp~^(y a ),P) E J ' u(x a ), (aexp~^(x a ),Q) 6 J ' v(y a ). Since u is 
subsolution and v is supersolution, and -F is continuous we then have, according to 
Remark l2~12l that 

F(xa,,u(x a ),—aexp~*(y a ),P) < < F(y a ,v(y a ),aexp~*(x a ),Q). (6) 

By combining equations (3), (4), (5) and (6) above, and using condition (1) too, 
we finally get 

< 7^ < 7 (u(x a ) - v(y a )) < 

F(x a ,u(x a ), -aexp~*(y a ), P) ~ F(x a , v(y a ), -aexp~*(y a ), P) < 
F(x a ,u(x a ), -aexp~^(y a ), P) - F(y a ,v(y a ),aexp~^(x a ),Q) + 
+F(y a ,v(y a ),aev^>yl(x a ),Q) - F(x a , v(y a ), -aexp~*(y a ), P) < 
< Lu(ad(x a ,y a ) 2 + d(x a ,y a )), 
and the contradiction follows by letting a —> oo. □ 

Remark 4.3. Observe that, since aexp~ 1 (a:) = L xy (— aexp~ 1 (y)), condition (2) 
of Theorem 14.21 can be replaced with a stronger but simpler assumption, namely 
that 

F(y, r, L xy C, Q) - F(x, r, P) < u (ad(x, y) 2 + d(x, y)) 
for all x,y £ O, r E M, P E T 2jS {M) x ,Q £ T 2>g (M) v , C £ TM* satisfying (*). 

Remark 4.4. If we want to compare two solutions u and v of F = and we know 
that these functions are bounded by some R > (e.g. when M is compact) then it 
is obvious from the above proof that it suffices to require that conditions (1) and 
(2) of Theorem 14.21 be satisfied for all r, s in the interval [-R, R\. 

Proposition 4.5. If M has nonnegative sectional curvature, then condition (*) 
implies that P < L yx {Q). 
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Proof. Let Ai,...,A n be the eigenvalues of the restriction of A a to the subspace 
V = {(v,L xy v) : v 6 TM X } of TM X x TM y . By Proposition 15711 we have that 
A a (v, L xy v) 2 < for all v G TM X , that is (A a )^ v < 0, or equivalently Xi < for 
i = 1, ...,n. With our choice of e a , this implies that 

A,: + e a X 2 < Xi + — - pr— tt-Aj <A, + ^ = ^<0 

2(l + aup 1 < J .< n |A J -|) 2 2 

and since Xi + e a Af , i = 1, n, are the eigenvalues of (A Q + e a A^) this means 
that 

v) 2 < 0. 

Then condition (*) implies that 

P{vf - Q{L xy vf < (A a + e a A 2 a ) {v, L xy vf < 
for all v £ TM X , which means that P < L yx (Q). □ 

Therefore, if M has nonnegative curvature and F is degenerate elliptic then (*) 
automatically implies that 

F(x,r,t,L yx Q)-F(x,r,(,P) < 0, 

hence 

F(y,r,L xy t,Q)-F(x,r,t,P) = 

F(y, r, L xy (, Q) - F(x, r, £ L y;c Q) + F(or, r, £ L^Q) - r, C, P) < 
P(?/, L^O Q) - r, C, L yx Q), 

and we see that condition (2) of the Theorem is satisfied if we additionally require, 
for instance, that 

F(y, r, rj, Q) ~ F(x, r, L yx n, L yx Q) < uj(d(x, y)). (2jt) 

Note that, in the case M = M™ we have L yx r) = rj and L yx Q = Q, and condition 
(2(t) simply means that F(x, u, du, d 2 u) is uniformly continuous with respect to the 
variable x. Therefore we can regard condition (2jJ) as the natural extension to the 
Riemannian setting of the Euclidean notion of uniform continuity of F with respect 
to x. This justifies the following 

Definition 4.6. We will say that F : X — * R is intrinsically uniformly continuous 
with respect to the variable x if condition (2fl) above is satisfied. 

Remark 4.7. As we saw in Example 12. f 41 above, many interesting examples of 
equations involving nonincreasing symmetric functions of the eigenvalues of d 2 u 
(such as the trace and the positive determinant det+) automatically satisfy condi- 
tion (2(() as long as they do not depend on x. In fact, since the eigenvalues of A 
are the same as those of L yx AL xy , any function of the form 

F(x, r, C,A) = G(r, ||C|Uj eigenvalues of A) 

is intrinsically uniformly continuous with respect to x. 

Therefore, for manifolds of nonnegative curvature, we do not need to impose 
that F depends on d 2 u{x) in a uniformly continuous manner: the assumptions that 
F is degenerate elliptic and intrinsically uniformly continuous with respect to x are 
sufficient. Let us sum up what we have just shown. 
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Corollary 4.8. Let £1 be a bounded open subset of a complete finite- dimensional 
Riemannian manifold M with nonnegative sectional curvature, and F : X — > K be 
continuous, degenerate elliptic, and satisfy: 

(1) F is strongly increasing, that is there exists 7 > such that, if r > s then 

1 (r-s)<F(x,r,t,Q)-F(x,s,<;,Q); 

(2) F is intrinsically uniformly continuous with respect to x (that is there exists 
a function lo : [0, 00] — > [0, 00] with lim t _, + u>(t) = and such that 

F(y, r, r,, Q) - F(x, r, L yx C, L yx Q) < u(d{x, y)) (2(t) 
for allx,y,r,(,Q with d(x,y) < mm{i M (x),i M (y)}). 
Let u G USC(Q) be a subsolution and v G LSC(Q) a supersolution of F = on fl, 
and u < v on d£l. 

Then u < v on all of fl. 

In particular, the Dirichlet problem (DP) has at most one viscosity solution. 

When M has negative curvature, condition (*) does not imply P < L yx Q, and 
degenerate ellipticity together with fulfillment of (2jJ) is not enough to ensure that 
condition (2) of Theorem 14.21 is satisfied. In this case condition (2) of 14. 21 involves 
kind of a uniform continuity assumption on the dependence of F with respect to 
d 2 u(x). Let us be more explicit. 

Proposition 4.9. Assume that M has sectional curvature bounded below by some 
constant —Kq < 0. Then condition (*) in Theorem \4-.2\ implies that 

P-L yx (Q) < ^K Q ad(x,y) 2 I, 

where I(v) 2 — (v,v) — \\v\\ 2 . 

Proof. We have that A a = (a/2)d 2 ip(x,y), where ip(x,y) = d(x,y) 2 . According to 
Proposition 13 . 31 we have 

d 2 <p(x,y)(v,L xy v) 2 < 2K d(x,y) 2 \\v\\ 2 

for all v G TM X and x,y £ M with d(x,y) < min{iM(^), %m (y)}- Therefore 

A a (v,L xy v) 2 < aK d(x,y) 2 \\v\\ 2 . 

This means that the maximum eigenvalue of the restriction of A a to T> :— {(v, L xy v) : 
v G TA1 X }, which we denote A n , satisfies 

A„ < aK d(x,y) 2 . 

If Ai, A„ are the eigenvalues of (A a )\D then + e a A?, i = 1, n, are those of 
(A a + Ea^jp. For a given i = 1, ...,n, if Aj < then A, + e a X 2 < as in the 
proof of of Remark 14.51 In particular, if A n < then Aj < for all j = 1, n, so 
we get that Aj + e a X 2 < for all j = 1, n, which means that 

(A a +e a A 2 a ) (v,L xy v) 2 < 0. 

On the other hand, if A„ > then A„ + £ a X 2 % > 0, and because the function 
[0, +00) 3 s 1— > s + e a s 2 G [0, +00) is increasing, the maximum eigenvalue of 
(A a + SaA^)^ is precisely A„ + e Q A^. This means that 

(A a + e a A 2 a ) (v, L xy v) 2 < A„ + e a \ 2 n 
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for all v G TM X . Besides we have, by the choice of e a , that 




hence 



(A a +e a A 2 a ) (v,L xy v) 2 < -A» < -aK d{x,y) 2 \\v\\ 2 . 



In any case (no matter what the sign of A„ is) we get that the above inequality 
holds. Therefore condition (*) implies 



Corollary 4.10. Let M be a complete Riemannian manifold (no assumption on 
curvature) , and SI be a bounded open subset of M . Suppose that F : X — > K is 
proper, continuous, and satisfies the following uniform continuity assumption: for 
every e > there exists 6 > such that 



d(x, y)<5, P- L yx (Q) <SI => F(y, r, L xy Q, Q) - F(x, r, ( 1 P)<£ (2b) 

for all x,y g M with d(x,y) < i n , r G R, ( € TM*, P e rf^j^^^ and 
Q £ C 2 (TM y ,M). Assume also that there is 7 > such that 



Then there is at most one viscosity solution of the Dirichlet problem (DP). 

Proof. Since M is complete, we know from the Hopf-Rinow Theorem that Q is 
compact, hence we have that in — inf^^iM^) > 0. Take a number r with 
< 2r < inf^gQ im(x). Also by compactness of f2, there exists K > such that 
the sectional curvature is bounded below by —Kq on f2. Therefore we have that 
<p(x,y) = d(x,y) 2 is C°° smooth on the set {(x,y) G ft x ft : d(x,y) < r} and, 
according to the preceding Remark, if P, Q satisfy condition (*) of Theorem 14.21 
we get P — L yx (Q) < ^Koad(x,y) 2 I whenever d(x,y) < r. Then the uniform 
continuity assumption on F yields the existence of a function to : [0, 00] — > [0, 00] 
with lim t ^ + w(t) = and such that 



F(y,r,aexp y 1 (x),Q) - F(x,r,-aexp x 1 (y),P) < u (ad(x,y) 2 +d(x,y)) , 



Remark 4.11. As is usual with comparison principles, the proof of Theorem 14.21 
can easily be adapted to show that the viscosity solutions u of the equations F = 
depend continuously on F . That is, if u is solution of F = and v is solution of 
G = 0, then 

sup \u(x) — v(x)\ < sup \F(x, r, (, A) — G(x, r, (, A)\ + sup \u — v\. 




□ 



j(r -s)< F(x, r, (, Q) - F(x, s, C, Q) for all r > s. 



hence the result follows from Theorem 14.21 



□ 
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5. Comparison results without boundary conditions 

The same argument as in the proof of Theorem 14. 2\ with some small changes, 
yields the following. 

Theorem 5.1. Let M be a connected, complete Riemannian manifold (without 
boundary) such that i(M) > 0, and F : X — > K be proper, continuous, and sat- 
isfy assumptions (1) and (2) of Theorem \4-2\ Let u be a subsolution, and v a 
super solution, of F — 0. Assume that u and v are uniformly continuous and 
linig;— too u(x) — v(x) < 0. Then u < v on M. In particular, if M is compact, 
there is at most one viscosity solution of F = on M. 

Uniform continuity and the inequality at infinity guarantee that the m a are 
attained, so the only difference with the proof of Theorem 14.21 is that now we 
cannot assume that x a and y a converge to some point xq, but we do have that 
d(x a ,y a ) < i(M) for large a, hence all the computations and estimations in the 
proof of 14. 21 are still valid. 

Corollary 5.2. Let M be a compact Riemannian manifold of nonnegative sectional 
curvature, and let F : X — > R be continuous, degenerate elliptic, strongly increasing 
in u, and intrinsically uniformly continuous with respect to x (that is, F satisfies 
conditions (1 — 2) or Corollary \4-8\ Let u be a subsolution, and v a supersolution 
of F = 0. Then u < v on M . 

Proof. The same considerations as in Remark 14 . 5 1 apply. □ 

Example 5.3. As we remarked above, condition (2) of Corollary 14.81 is easily 
satisfied when F(x, r, £, A) does not depend on ( and A themselves, but on ||£|| and 
the eigenvalues of A. For instance, the function 

F{x, r, C, A) = r-(dct + (A)) 3 ||C|| 2 - /(z)(trace(A)) 5 

satisfies (1) and (2) of the above Corollary provided that / > and / is uniformly 
continuous. Therefore the equation 

u - (Act+{D 2 u)f || Vu|| 2 - (Au) 5 f = 

has at most one viscosity solution on any compact manifold of positive curvature 
if we only require that / is continuous and nonnegative. 

Corollary 5.4. Let M be a compact Riemannian manifold (no assumption on 
curvature). Suppose that F : X — > M satisfies the uniform continuity assumption, 
and the growth assumption, of Corollary \4-10\ Let u be a subsolution, and v a 
supersolution of F = 0. Then u < v on M . In particular there is at most one 
viscosity solution of F = 0. 

Proof. Since M is compact the sectional curvature of M is bounded on all M, say 
K > —Kq. Take a number r with < 2r < i(M). The function tp(x,y) — d(x,y) 2 
is C°° on the set {(x, y) £ M x M : d(x, y) < 2r}. Suppose that P and Q satisfy (*) 
of Theorem |4*T21 Then, from Remark WM we get that P - L yx Q < \Kq a d(x, y) 2 I 
provided that d{x, y) < r. Therefore the uniform continuity property of F gives us 
a function lj : [0, oo] — > [0, oo] with lim t ^ + ^(t) =0 and such that 

F{y,r,aexp~ l {x),Q) - F(x,r,-aex$~ x (y),P) <u(ad(x,y) 2 +d(x,y)) . 

Hence we can apply Theorem 15.11 and conclude the result. □ 
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6. Existence results 

Perron's method can easily be adapted to the Riemannian setting to establish 
existence of viscosity solutions to the Dirichlet problem. The proof goes exactly as 
in [4] with appropriate changes. The only step which is not completely obvious is 
the proof of the following 

Proposition 6.1. Let (£, A) G J 2 ' + f(z) Suppose that f n is a sequence of upper 
semicontinuous functions such that 

(i) there exists x n such that {x ni f n (x n )) — > (x,f(x)), and 

(ii) ifVn ->■ V, then limsup^^ f n (y n ) < f(y). 

Then there exist x n and (£ n ,A n ) G J 2,+ f n (x n ) such that (x n f n (x n ),( n ,A n ) — » 
(xJ(x)X,A). 

Proof. Consider the functions / o exp^, and /„ o exp^, defined on a neighborhood of 
in TM X . These functions satisfy properties (i) and (ii) of the statement (when 
they take the roles of / and /„ and M is replaced with TM X ). By Corollary 12. 31 we 
have that ((, A) G J 2,+ (f ° exp a ,)(0). And of course the result is known in the case 
when M = W 1 , so we get a sequence v n and (( n ,A n ) G J 2 ' + (f n ° exp a .)(u„) such 
that 

(^nj fn 

oexp x (v n ),(n,A n ) -> (0,foexp x (0),(,A). 

Set x n = exp x (v n ). We have that x n — > x and f n (x n ) — * f{x). Since (( n ,A n ) G 
J 2 ' + {fn ° exp^)^) there exist functions ip n such that /„ o exp x —ip n attains a 
maximum at v n , ( n = dip n (y n ) and A n = d 2 ip„(v n ). Let us define ip = ip n o exp" 1 
on a neighborhood of x. Then /„ — (p n attains a maximum at x n so, if we set 
C n = dip(x n ), A n = d 2 ip(x n ), we have that A n ) G J 2,+ f n {x n )- It only remains 
to show that £ n — > ^ and A n — > A. But this is exactly what was shown in (<=) of 
the proof of Proposition 12.81 □ 

By using this Proposition one can prove, as in [4] , existence of viscosity solutions 
to the Dirichlet problem 

F(x, u, du, d 2 u) = in ft, u = f on <9ft, (DP) 

where 51 is an open bounded subset of a complete Riemannian manifold M. 

Theorem 6.2. Let comparison hold for (DP), i.e., if w is a subsolution of (DP) 
and v is a supersolution of (DP), then w < v. Suppose also that there exists a 
subsolution u and a supersolution u of (DP) that satisfy the boundary condition 
u*(x) — u*(x) — f(x) for x G dfl. Then 

W(x) — sup{w(x) : u< w <u and w is a subsolution of (DP)} 

is a solution of (DP). 

Here we used the following notation: 

u*(x) = lim r ^o sup{u(y) : y G fl and d(y,x) < r}; 
u*(x) — lim r jo inf {u(y) ■ y E and d(y,x) < r}, 

that is u* denotes the upper semicontinuous envelope of u (the smallest upper 
semicontinuous function, with values in [—00,00], satisfying u < u*), and similarly 
it* stands for the lower semicontinuous envelope of u. 
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One can also easily adapt the proof of [21 Theorem 6.17] to the second order 
situation, obtaining the following. 

Corollary 6.3. Let M be a compact Riemannian manifold, and G(x, du, d 2 u) be 
degenerate elliptic and uniformly continuous in the sense of Corollary \^.1C\ Then 
there exists a unique viscosity solution of the equation u + G(x,du,d 2 u) — on 
M. 

Again, if M has nonnegative curvature, the assumptions that F is elliptic and 
intrinsically uniformly continuous with respect to x are sufficient in order to get an 
analogous result. 

7. Examples 

Most of the examples of proper F's given in [4] remain valid in the Riemann- 
ian setting. In particular, as we have already seen, the functions (x, r, £, A) >— > 
— det + (A) and (x, r, £, A) t— > — trace(A) are degenerate elliptic and intrinsically uni- 
formly continuous with respect to x. The same is true of all many symmetric 
functions of the eigenvalues of A, such as minus the minimum (or the maximum) 
eigenvalue, and of course nondecreasing combinations and sums of these are degen- 
erate elliptic too. One can find lots of examples of nonlinear equations for which the 
results of this paper yield existence and uniqueness of viscosity solutions. For in- 
stance, one can easily show that, for every compact manifold of positive curvature, 
the equation 

max{u - Ai(D 2 u)||Vu|| p - (Am) 2 « +1 || Vuf - (det + (£> 2 it)) 2fc+1 f , u - g} = 

(where Ai denotes the minimum eigenvalue function and p, q, r, k G N) has a unique 
viscosity solution if we only require that / and g are continuous. This gives an idea 
of the generality of the above results. 

Of course this example is rather unnatural. Let us finish this paper by examining 
what our results yield in the case of a classic equation, that of Yamabe's, which 
has been extensively studied and completely solved by using variational methods. 
We do not claim that the following discussion gives any new insight into Yamabe's 
problem, we only want to study, from the point of view of the viscosity solutions 
theory, a well known example of a nonlinear equation arising from an important 
geometrical problem. 

Example 7.1 (The Yamabe equation). A fundamental problem in conformal ge- 
ometry is to know whether or not there exists a conformal metric g' with constant 
scalar curvature S' on a given compact n-dimensional Riemannian manifold (M, g), 
with n > 3, see [II I12j. This is equivalent to solving the equation 

77 — 1 n-t-2 

-4- Au + S(x)u = S'u— , (Y) 

rt — 2 

where S is the scalar curvature of g. One can write this equation in the form F = 0, 
where 

n — 1 

Fix, r, C, -4) = S(x)r - 5'r~ - 4 trace(A) = 0. 

n — 2 

It is clear that F is degenerate elliptic. Assume that S is everywhere positive and 
that S' < 0. Then, by compactness, there exists 7 > such that S(x) > 7 for 
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all x € M. According to Remark 14. 4\ in order to check conditions (1) and (2) of 
Theorem 15.11 we may assume that r, s lie on a bounded interval. We have that 

F(y,r,r),Q) - F(x,r,Ly X i],L yx Q) < r\S(y) - S(x)\, 

hence, because S is uniformly continuous on M and r is bounded, we deduce that 
F satisfies (2) of Corollary [521 On the other hand, if r > s then 

F(x,r,(,A) -F(x,s,C,A) = S(x)(r - s) - S'(r^ - s^) > j(r-s), 

so condition (1) is also satisfied. It follows that there is at most one viscosity 
solution of F = 0. Existence can be shown by using Perron's method. In all, we 
see that if 5* is everywhere positive and S' < then there exists a unique viscosity 
solution u of (Y). 
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